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A Roe-type Riemann solver based on the spectral decomposition of 
the equations of Relativistic Magnetohydrodynamics 



Jose M- Ibanez' , Miguel A. Aloy \ Petar Mimica' , Luis Anton' , 
Juan A. Miralles^, Jose Marti' 

Abstract. In a recent paper dAnton et al.ll2010h we have derived sets of right and left 
eigenvectors of the Jacobians of the relativistic MHD equations, which are regular and 
span a complete basis in any physical state including degenerate ones. We present a 
summary of the main steps followed in the above derivation and the numerical exper- 
iments carried out with the linearized (Roe-type) Riemann solver we have developed, 
and some note on the (non-)convex character of the relativistic MHD equations. 



1. Introduction 

Relativistic flows in association with intense gravitational and magnetic fields are com- 
monly linked up to extremely energetic phenomena in the Universe, viz. pulsar winds, 
anomalous X-ray pulsars, soft gamma-ray repeaters, gamma-ray bursts, relativistic jets 
in active galactic nuclei, etc. The necessity to model the aforementioned astrophysical 
scenarios in the framework of relativistic MHD (RMHD), together with the fast increase 
in computing power, is pushing towards the development of more efficient numerical 
algorithms. In the last years, considerable progress has been achieved in numerical 
special RMHD (SRMHD), by extending the existing high-re solution shock- capturing 
(HRSC) methods of special relativistic hydrodynamics (e.g.. lMartf & Mulleril2003b . In 
the so called Godunov-type methods, an important subsample of HRSC methods, nu- 
merical fluxes are evaluated through the exact or approximate solution of the (local) 
Riemann problem. Despite the fact that such an ex act solution in SRMHD is known 
( Romero et al.ll2005l : iGiacomazzo & RezzoUal 12006 '). approximate algorithms are usu- 



ally preferred bec ause of their larger numerical efficiency. Several authors (see, e.g., 
Anton et aLll2010l and references therein) have developed independent Roe-type algo- 



rithms based on linearized Riemann solvers relying on the characteristic structure of 
the RMHD equations. 

The purpose of the present paper is twofold. On one hand, the objective is to 
present a regular set of right and left eigenvectors of the flux vector Jacobian matrices 
of the RMHD equations, and span a complete basis in any physical state, including 
degenerate states. On the other hand, wish to evaluate numerically the performance of 
a RMHD Riemann solver based on the aforementioned spectral decomposition. Both 
the theoretical analysis and the numerical applications presented in this paper are based 



'Department of Astronomy and Astrophysics, University of Valencia, 46100 Burjassot (Valencia), Spain 
"Department of Applied Physics, University of Alicante, Ap. Correus 99, 03080 Alacant, Spain 

I 



2 



on the work developed by lAnton et alj ( 20Toh . where we have characterized thoroughly 
all the degeneracies of RMHD in terms of the components of the magnetic field nor- 
mal and tangential to the wavefront in the fluid rest frame. Our numerical method 
deviates in several aspects fror t i previoiis wo r ks based on linearize d Riemann solver 
approaches ('Komissarovl ll999l : lBalsara'l 200ll : iKoldoba et"aill2002l ). First, numerical 
fluxes are computed from the spectral decomposition in conserved variables. Second, 
we present explicit expressions also for the left eigenvector s. Third, and most impor- 
tant, we have extended classical MHD strategy ('Brio & Wu"l988) to relativistic flows, 
giving sets of right and left eigenvectors which are well defined through degenerate 
states. Based on the full wave decomposition (FWD) provided by the renormalized set 
of eigenvectors in conserved variables, we have also developed a linearized (Roe-type) 
Riemann solver. 

Extensive testing against one- and two-dimensional standard numerical problems 
allows us to conclude that our solver is very robust. When compared with a family of 
simpler solvers that do not require the knowledge of the full characteristic structure of 
the equations in the computation of the numerical fluxes, our solver turns out to be less 
diffusive than HLL and HLLC, and comparable in accuracy to the HLLD solver. The 
amount of operations needed by the FWD solver makes it less efficient computation- 
ally than those of the HLL family in one-dimensional problems. However its relative 
efficiency increases in multidimensional simulations. 



2. The equations of ideal relativistic magnetohydrodynamics 

The equations of ideal RMHD correspond to the conservation of rest-mass and energy- 
momentum, and the Maxwell equations. In the following, the standard Einstein sum 
convention is assumed. Greek indices will run from to 3 (or from t to z) while Roman 
run from 1 to 3 (or from x to z). We use units in which the speed of light is c = 1 and 
the (47r)^^^ factor is absorbed in the definition of the magnetic field. Specializing for a 
flat space-time and Cartesian coordinates, these equations can be written as a system of 
conservation laws, which reads 

— + — r - 0, (1) 

dt dx' 

where the state vector, U (vector of conserved variables), and the fluxes, F' (/ = 1, 2, 3 
or / = x,y,z), are the following column vectors, 

U - (D,5^',r,By (2) 

F' = {Dv\ S + p*5'^ - yB'IW, W + p*V - b^B'/W, Vb'' - v'^B'f (3) 

where the superscript ^ stands for the transposition. 

In the preceding equations, D, and r stand, respectively, for the rest-mass den- 
sity, the momentum density of the magnetized fluid in the j-direction and its total en- 
ergy density as measured in the laboratory (i.e., Eulerian) frame, 

D=pW, =ph*W^v} -b°bj , T = ph*W'^ - p* - {by -D. (4) 

where p is the proper rest-mass density, h* - l + e + p/p + b^/pis the specific enthalpy 
including the contribution from the magnetic field (b^ stands for b^b^), e is the specific 
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internal energy, p the thermal pressure, and p* = p + /2 the total pressure. The four- 
vectors representing the fluid velocity, m'^, and the magnetic field measured in the fluid 
rest frame, b^^, and there is an equation of state relating the thermodynamic variables, 
p, p and e, p = p{p, e). All the discussion will be valid for a general equation of 
state but results will be shown for an ideal gas, for which p = {y - l)pe, where y is the 
adiabatic exponent. Quantities v' stand for the components of the fluid velocity trivector 
as measured in the laboratory frame; they are related with the components of the fluid 
four-velocity according to the following expression u^^ = W(l,v^, v-^, v^), where W is 
the flow Lorentz factor, = 1/(1 - v'v,). 

The following fundamental relations hold between the components of the magnetic 
field four-vector in the comoving frame and the three vector components B' measured 
in the laboratory frame, 

b^^WB-\ , b' = — + b^V (5) 
W 

V and B being, respectively, the tri-vectors (v^, v-^', v^) and (B^, B^, B'). 

The preceding system must be complemented with the usual divergence constraint 

dB' 

which should be fulfilled at all times. 

Fluxes F' (/ = X, y, z) are functions of the conserved variables, U, although for the 
RMHD this dependence, in general, can not be expressed explicitly. It is therefore nec- 
essary to introduce another set of variables, the so-called primitive variables, derived 
from the conserved ones, in terms of which the fluxes can be computed explicitly. We 
have used the following set of primitive variables 

\ = (p,p, v\ v^, v^ B\ By, B'f. (8) 



3. Characteristic structure of the RMHD equations 

The hyperbolicity of the equations of RMHD including the derivation of wavespeeds 
and the corr esponding e igenvectors, and the analysis of various degeneracies has been 
reviewed by lAnild (Il989 ). in a covariant framework, using a set of variables of dimen- 
sion 10, the so-called covariant variables (Anile's variables, in the next): 

I] = {u'',b'',p,sf, (9) 

where s is the specific entropy. 

In terms of variables U, the system of RMHD equations can be written as a quasi- 
linear system of the form 



^''U;/, - 0, 



(10) 
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where the subscript ;/i stands for the covariant derivative, and four 10 x 10 Jacobian 
matrices can be found in Anile's book. It is important to remark that the 10 covariant 
variables we have used to write the system of equations are not independent, since they 
are related by the constraints 

u"ua = -1 , b^Ua = , daiu^b^ - u^b") - 0, (11) 
The latter condition, is a covariant representation of the divergence constraint (Eq.|7]l. 

3.1. Wavespeeds and degeneracies 

The system of (ideal) RMHD equations have the same seven wavespeeds as in classical 
MHD: the entropic, Alfven, slow magnetosonic, and fast magnetosonic waves. They 
can be ordered as follows 

A} < A- <A-<Ae<A:< Al < (12) 

where the subscripts e, a, s and / stand for entropic, Alfven, slow magnetosonic mAfast 
magnetosonic respectively, and the superscript - or + refer to the lower or higher value 
of each pair. Unlike classical MHD, it is however not possible, in general, to obtain 
simple expressions for the magnetosonic speeds since they are given by the solutions 
of a quartic equation. 

As in the case of classical MHD, degeneracies are encountered for waves propa- 
gating perpendicular to the magnetic field direction (Type I) and for waves propagating 
along the magnetic field direction (Type II). Finally, a particular subcase of Type II 
degeneracy appears when the sound speed is equal to Ca- 

For Type I degeneracy, the two Alfven waves, the entropic wave and the two slow 
magnetosonic waves propagate at the same speed (i^ = A''^ = Ag = A'^ = A'^). For Type 
II degeneracy, an Alfven wave and a magnetosonic wave (slow or fast) of the same 
class propagate at the same speed {A~j = A~ or A~ = A~ or A'^ - A'^ or A'^ = ip. In 
the special Type IF subcase, an Alfven wave and both the slow and fast magnetosonic 
waves of the same class propagate at the same speed (AJ = A~ = A~ or A^ = A'^ = A'^). 

3.2. Renormalized right eigenvectors 

As it is well known in classical MHD, Alfven and magnetosonic eigenvectors have a 
pathological behaviour at dege neracies, since they become zero or linearly dependent 
and they do not form a basis. In lAnton et al. (I2010h . we have derived a new set of renor- 



malized Alfven and magnetosonic eigenvectors for RMHD. Our renormalized Alfven 
right eigenvectors (following Bri o & Wa,198 8 methodology) are a linear combination 
of the ones proposed by Komissarovl ( 1999h . for the Type II degeneracy case. However, 



contrary to the Komissarov's choice, our expressions are free of pathologies not only 
in the Type II degeneracy but also in the Type I degeneracy case. Our derivation of 
the renormalized magnetosonic righ t eigenvectors is al gebraically more cumbersome 
and reader interested is addressed to lAnton et"aL I teoioh . The final result of this anal- 



ysis allows to have a complete set of right eigenvectors linearly independent for all 
possible states. Following the sam e procedure we ha ve used to renormalize the right 
eigenvectors, we have derived (see lAnton et al.ll2010l) left eigenvectors well behaved 
for degenerate states. 
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4. A Full Wave Decomposition Riemann Solver in RMHD (FWD) 

Let us summarize the main steps allowing us to derive a FWD Riemann Solver in 
RMHD: 

• Let r~ be a generic right eigenvector derived in terms of Anile 's variables 

• The corresponding eigenvector in terms of the primitive variables ^ is derived 
according to: = (d^\^r^. 

• Finally, the corresponding vector in terms of the conserved variables Q is ob- 
tained from R = = id\U) r^. 

• Analogously, for the left eigenvectors. This procedure, which starts with renor- 
malized eigenvectors, allows one to get the full spectral decomposition of the 
Jacobian matrices (associated to the fluxes), and free of pathologies in the degen- 
eracies. 

• We use a hnearized (Roe's type) Riemann solver: 

f/±i - h (f("t 1 ) + f(ul 1 ) - K-i I 1 ^^«^-] 

where u^, u^, are the left and right reconstructed variables; Aw, is the jump of 
characteristic variables. 



^Our FWD linearized Riemann solver has been exhaustively tested (lAnton et al. 

2010h . 



5. A note on RMHD convexity 

For the sake of conciseness let us remind some definitions. A characteristic field 
Ca{a - 1, 2, iif) (cf is the number of equations) satisfying 

Ca- ^ = ^a {a=\,2,...,d) (13) 
at 

is said to be genuinely nonlinear or linearly degenerate if, respectively, 

Vu ^ • r„ ^ 0, (14) 
Vu ^ • r„ - (15) 

where the operator Vu acts on the space of conserved variables. 

In a convex system, all the characteristic fields are genuinely non-linear or linearly 
degenerate. Non-convexity is associated to those states for which the condition ([T4l ) is 
not fulfiled. 

For the system of equations governing relativistic (ideal) flows it can be shown that 
the convexity is strongly dependent on the second derivatives of pressure (or the first 
derivatives of the sound speed). In the following, we analize when genuinely nonlinear 

fields become linearly degenerate, by examining the products f + := V„ A±{v/) ■ r+(w). 
After some algebra, we find 

^,= ±r(«,Z,,c.)(^ + 4^ + ££(l_c2)) (16) 
\op p-^ oe p I 
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being 

T{a,b,c,,) = (ac, ± S^'^T^ (l - a^f d~^'^ (17) 

where a and b stand for, respectively, the spatial components of the velocity field in the 
x-direction and the tangential one. The quantity 6 is defined hy 6 = W^{1 - - b^c]). 

^From the above relations ( I16I17I ) it turns out that the loss of convexity is closely 
related with the properties of the equation of state (second term in Eq.[T6l). The first and 
second thermodinamical deriv atives of pressure p l ay a fu ndamental role regarding with 
this issue (that was noticed by iMenikoff & Pldij (Il989h . for equations of state having 
phase transitions). Furthermore, we realize that 

i) In the purely one-dimensional case {b = 0), non-convexity only appears in the ultra- 
relativistic regime (a — > 1): 

T{a,0,c,) = {\ ± ac-y^W--^ ■ a ^ \ ^ T{a,0,c,) 0, 0(W"^) (18) 

ii) Likewise, if a = 0, non-convexity arises in the ultrarelativistic regime (b — > 1): 

T{0,b,Cs) - A~^'^ ; b ^ I ^ T{0,b,c,) 0, OiW'^) (19) 
Brio & 1988 ) noted that the equations of classical MHD are non-convex at the 



degenerate states (magnetosonic waves change from genuinely non-linear to linearly 
degenerate) . We have faced on the problem of non-convex character of RMHD, and 
preliminar results allows one to conclude that the degenerate states are, as in the classi- 
cal MHD, non-convex , being t he ma gnetosonic fields the ones changing their character. 
We refer the reader to lAntonI ( 20081) (appendix G), where an analysis of the character- 



istic fields of RMHD in terms of Anile's covariant variables is presented. Much more 
theoretical work is necessary in order to asses all the richness of other possible non- 
convex states in RMHD. The previous analysis in special relativistic hydrodynamics 
serves us as a road-map to the full characterization of RMHD pathological behaviours 
(we remind the reader that the non-convex character of both the classical and relativistic 
MHD equations is source of several pathologies, as the development of the so-called 
compound waves). 
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